Dynamics of collapsing and exploding Bose-Einstein condensates 
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We explored the dynamics of how a Bose-Einstein condensate collapses and subsequently explodes 
when the balance of forces governing the size and shape of the condensate is suddenly altered. A 
condensate's equilibrium size and shape is strongly affected by the inter-atomic interactions. Our 
ability to induce a collapse by switching the interactions from repulsive to attractive by tuning 
an externally-applied magnetic field yields a wealth of detailed information on the violent collapse 
process. We observe anisotropic atom bursts that explode from the condensate, atoms leaving the 
condensate in undetected forms, spikes appearing in the condensate wave function, and oscillating 
remnant condensates that survive the collapse. These all have curious dependencies on time, the 
strength of the interaction, and the number of condensate atoms. Although ours would seem to 
be a simple well-characterized system, our measurements reveal many interesting phenomena that 
challenge theoretical models. 



Although the density of the atoms in an atomic Bose- 
Einstein condensate (BEC) is typically five orders of mag- 
nitude lower than the density of air, the inter-atomic in- 
teractions greatly affect a wide variety of BEC properties. 
These include static properties like the BEC size and 
shape and the condensate stability, and dynamic prop- 
erties like the collective excitation spectrum and soliton 
and vortex behavior. Since all of these properties are sen- 
sitive to the inter-atomic interactions, they can be quite 
dramatically affected by tuning the interaction strength 
and sign. 

The vast majority of BEC physics is well described by 
mean-field thcoryu, in which the strength of the interac- 
tions depends on the atom density and on one additional 
parameter called the s-wave scattering length a. a is 
determined by the atomic species. When a > 0, the in- 
teractions are repulsive. In contrast, when a < the 
interactions are attractive and a BEC tends to contract 
to minimize its overall energy. In a harmonic trap, the 
contraction competes with the kinetic zero-point energy, 
which tends to spread out the condensate. For a strong 
enough attractive interaction, there is not enough kinetic 
energy to stabilize the BEC and it is expected to implode. 
A BEC can avoid implosion only as long as the .number 
of atoms No is less than a critical value given bya 



N cr = kaho/\a\ , 



(1) 



where dimensionless constant k is called the stability co- 
efficient. The precise value of k depends on the aspect 
ratio of the magnetic trapcl. aho is the harmonic oscilla- 
tor length, which sets the size of the condensate in the 
ideal-gas (a = 0) limit. 

Under most circumstances, a is insensitive to external 
fields. This is different in the vicinity of a so-called Fesh- 
bach resonance, where a can be tuned over a huge rar. 
by adjusting the externally applied magnetic fieldE 
This has been demejftsfcrated in recent years with cold 
85 Rb and Cs atomiffl and with Na and 85 Rb Bose- 
Einstein condensatesBEI For 85 Rb atoms, a is usually 
negative, but a Feshbach resonance at ~155 G allows 



us to tune a by orders of magnitude and even change 
its sign. This gives us the-ability to create stable 85 Rb 
Bose-Einstein condensates^ and adjust the inter-atomic 
interactions. We recently used this flexibility to verify 
the functional form of equation (Q) pad to measure the 
stability coefficient to be k = 0.46(6)0. 

In this article, we study the dynamical response ("the 
collapse") of an initially stable BEC to a sudden shift 
of the scattering length to a value more negative than 
the critical value a cr = —kai lo /No- We have observed 
many features of the surprisingly complex collapse pro- 
cess, including the energies and energy anisotropies of 
atoms that burst from the condensate, the time scale for 
the onset of this burst, the rates for losing atoms, spikes 
in the wave function that form during collapse, and the 
size of the remnant BEC that survives the collapse. The 
unprecedented level of control provided by tuning a has 
allowed us to see how all of these quantities depend on 
the magnitude of a, the initial number and density of 
condensate atoms, and the initial spatial size and shape 
of the BEC before the transition to instability. 

A great deal of theoretical interestEiS330O was 
generated by BEC experiments in 7 LiE£l, for which the 
scattering length-is also negative and collapse events are 
also observedt£lEj. The 7 Li experiments do not employ 
a Feshbach resonance, so a is fixed. This restricts ex- 
perimentation to the regime where the initial number of 
condensate atoms is less than or equal to N cr , and the 
collapse is driven by a stochastic process. In addition, 
studies of collapse dynamics in 7 Li are complicated by a 
large thermal component. Our ability to tune the scatter- 
ing length, and most notably explore the regime where 
the initial condensate is "pure" (near T = 0) and the 
number Nq is much larger than N cr , allows us to explore 
the dynamics and compare it with theory in far more 
detail. 

Experimental techniques. The procedure for pro- 
ducing stable 8 ^E,b condensates has been described in 
detail elsewhereEI -4, standard double magneto-optical 
trap (MOT) systemEII was used to collect a cold sample 
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of 85 Rb atoms in a low-pressure chamber. Once suffi- 
cient atoms had accumulated in the low-pressure MOT, 
the atoms were loaded into a cylindrically symmetric 
cigar-shaped magnetic trap with trap the frequencies 
Vradiai = 17.5 Hz and v axia i = 6.8 Hz. Radio-frequency 
evaporation was then used to cool the sample to ~3 nK 
to form pure condensates containing >90% of the sample 
atoms. The final stages of evaporation were performed at 
162 G where the scattering length is positive and stable 
condensates of up to 15,000 atoms could be formed. After 
evaporative cooling, the magnetic field was ramped adi- 
abatically to 166 G (except where noted), where a = 0. 
This provided a well-defined initial condition with the 
BEC taking on the size and shape of the harmonic oscil- 
lator ground state. 

We could then adjust the mean-field interactions 
within the BEC to a variety of values on time scales 
as short as 0.1 ms. The obvious manipulation was to 
jump to some value of a < a cr to trigger a collapse, but 
the tunability of a also greatly aided in imaging the sam- 
ple. Usually the condensate size was below the resolution 
limit of our imaging system (7/im FWHM). However, we 
could ramp the scattering length to large positive values 
and use the repulsive inter-atomic interactions to expand 
the BEC before imaging, thus obtaining information on 
the pre-expansion condensate shape and number. A typ- 
ical a(t) sequence is shown in Fig. |l|a. We have used a 
variety of such sequences to explore many aspects of the 
collapse and enhance the visibility of particular compo- 
nents of the sample. 

Condensate contraction and atom loss. When 
the scattering length is jumped to a value a co u apse < 
a cr , a condensate's kinetic energy no longer provides a 
sufficient barrier against collapse. As described in ref. [l^, 
during collapse one might expect a BEC to contract until 
losses from density-dependent inelastic collisionsEj would 
effectively stop the contraction. This contraction would 
roughly take place on the time scale of a trap oscillation, 
and the density would sharply increase after T rac [/4: ~ 
14 ms, where T ra d is the radial trap period. How does 
this picture compare to what we have actually seen? 

A plot of the condensate number N vs T evo ive for 
a-coiiapse = -30 a and a init = +7 a is presented in Fig. 
@b. N was constant for some time after the jump until 
atom loss suddenly began at t co u apse . After the jump 
the condensate was smaller than our resolution limit, 
so we could not observe the contraction directly. But 
we observed that the post-expansion condensate widths 
changed very little with time T evo i ve before t co u apse . From 
this we infer that the bulk BEC did not contract dramat- 
ically before loss began. We modelled the contraction by 
putting the a(t) sequence that we used for the experi- 
ment into the equations of ref. for the approximate 
evolution of the BEC shape. The model predicted that 
by tcoiiaps,,, the radial and axial widths had contracted 
by ^20% and ~4%, respectively, from their initial values. 
This contraction only corresponds to a 50% increase in 
the average density to 2.5 x 10 13 /cm 3 . Using the decay 



constants from ref. g2| this density gives an atom loss 
rate, Tdecay, that is far smaller than what we observe and 
does not have the observed sudden onset. 

For the data in Fig. [l]b and most other data presented 
below, we jumped to a qU ench = in 0.1 ms after a time 
Tevoive at a co u a pse- We believe that the loss immediately 
stopped after the jump. This interpretation is based on 
the surprising observation that the quantitative details 
of curves such as that shown in Fig. ijb did not depend 
on whether the collapse was terminated by a jump to 



or a 



250 oq. 



We have measured loss curves like that in Fig. [T}d 
for many different values of a co iiapse- The collapse time 
versus a co u apse for N = 6000 presented in Fig. g shows 
the strong dependence on a co iiapse- Reducing the initial 
density by a factor of ~ 4 (with a corresponding increase 
in volume) by setting = +89 ao for one value of 

^collapse 

(-15 a ), increased t coUapse by — 3. 
The atom loss time constant Tdecay depended only 
weakly on a co Uavse and A . For the range of a co u apse 
shown in Fig. g, T deC ay did not depend on a coUapse or 
Ao outside of the experimental noise (~20%). On aver- 
age, Tdecay was 2.8(1) ms. For the very negative value 



of 



^collapse 



-250 ao, however, Tdecay did decrease to 



1.8 ms for A = 6, 000 and 1.2 ms for N = 15, 000. 

Burst atoms. As indicated by Fig. [j]b, atoms leave 
the BEC during the collapse. There are at least two 
components to the expelled atoms. One component (the 
"missing atoms") is not detected. The other component 
emerges as a burst of detectable spin-polarized atoms 
with energies much greater than the initial condensate's 
energy but much less than the magnetic trap depth. The 
burst-energy dependencies on a co iiapse and Ao are com- 
plex, but since they will provide a stringent test of col- 
lapse theories, we present them in detail. 

The angular kinetic-energy distribution with which the 
burst atoms are expelled from the condensate can most 
accurately be measured by observing their harmonic os- 
cillations in the trap, as illustrated in Fig. |3^,. For exam- 
ple, half of a radial period after the expulsion (T ra d/2), 
all atoms return to their initial radial positions. Well 
before or well after this "radial focus" , the burst cloud 
is too dilute to be observed. Fortunately, at the radial 
focus, oscillations along the axial trap axis are near their 
outer turning points and the axial energy can be found 
from the length of the stripe of atoms along the axial 
axis. The radial energy can be found with the same pro- 
cedure for an axial focus. The sharpness of the focus also 
provides information on the time extent of the burst. 

Figure [|b shows an image of a radial focus. The size 
scales for the burst focus and the remnant were well sep- 
arated since the latter was not expanded before imaging. 
Figure ||c shows cross sections of the burst focus and fits 
to the burst and the thermal cloud. The burst energy 
distributions were well fit by Gaussians characterized by 
a temperature that was usually different for the two trap 
directions. The burst energy fluctuated from shot to shot 
by up to a factor of 2 for a given a co iiapse- This varia- 
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tion is far larger than the measurement uncertainty or 
the variation in initial number (both ~ 10%), and its 
source is unknown. (We also discuss observed structural 
variability when we present the jet measurements below.) 

Although the burst energies varied from shot to shot, 
the average value was well-defined and showed trends far 
larger than the variation. The axial and radial burst 
energies versus a co ii a pse are shown in Fig. ^a and ^b for 
-/V = 6, 000 and 15,000, respectively. The burst-energy 
anisotropy shown in Fig. ^c depended on No, a co u apse , 
and a-init ■ 

When we interrupted the collapse with a jump back 
to a = as discussed above, we also interrupted the 
growth of the burst. The "interrupted" burst atoms still 
refocused after sitting at a — for the requisite half trap 
period. The energy of the atoms in the interrupted bursts 
appeared to be the same, but the number of atoms was 
smaller. By changing the time at which the collapse was 
interrupted we could measure the time dependence of the 
creation of burst atoms. For the conditions of Fig. |l|b, 
the number of burst atoms Nburst grew with T evo ive with 
a time constant of 1.2 ms starting at 3.5 ms and reaching 
an asymptotic final number of —2500 for all times >7 ms. 
Nburst varied randomly by —20% for the data in Fig. ||, 
but on average the fraction of atoms going into the burst 
was about 20% of iVo and did not depend on a co u apse or 
N . 

Remnant Condensate. After a collapse, a "rem- 
nant" condensate containing a fraction of the atoms sur- 
vived with nearly constant number for more than 1 sec- 
ond and oscillated in a highly excited collective state with 
the two lowest modes with v — 2v ax i a i and v — 2v ra( n a i 
being predominantly excited. (The measured frequencies 
were v = 13.6(6) Hz and v = 33.4(3) Hz.) To find the 
oscillation frequencies, the widths of the condensate were 
measured as a function of time spent at a co iiapse- 

The number of atoms in the remnant depended on 
o-coiiapse and No, and in general was not limited by the 
critical number, N cr . The stability condition in equa- 
tion (P determined the collapse point but it did not con- 
strain N remnant . A fixed fraction of No went into the 
remnant independent of iVo, so that smaller condensates 
often ended up with N remnant < N cr , but larger conden- 
sates rarely did. The fraction of atoms that went into 
the remnant decreased with |a co /; apse |, and was —40% 
for |a co/;apse | < 10 a and —50% for |oco/iop*e| > 1PQ «o- 
We do not think that the surface-wave excitationsc3 are 
responsible for stabilizing the remnant since we excite 
large-amplitude breathing modes. For No — 6, 000 and 
|flcoZ2apse| < 10 ao, more atoms were lost than the number 
required to lower N remnant to below N cr . 

Since N burst was independent of a coUapse but N remnant 
decreased with |a co /; apse |, the number of missing atoms 
increased with |a co « apse |. Interestingly, the number of 
missing atoms also increased with No , but the fraction of 
missing atoms versus acoiiapse was equal for No = 6, 000 
and N = 15,000 and was —40% for \a co u apse \ < 10 a 
and —70% for \a co u apse \ > 100 ao- The missing atoms 



were presumably either expelled from the condensate 
at such high energies that we could not detect them 
(> 20 /UK), or they were transferred to untrapped atomic 
states or undetected molecules. 

Jet formation. Under very specific experimental 
conditions, we observed streams of atoms with highly 
anisotropic velocities emerging from the collapsing con- 
densates. These "jets" are distinguished from the "burst" 
in that the jets have much lower kinetic energy (on 
the order of a few nanokelvin), in that their velocity 
is nearly purely radial, and in that they appear only 
when the collapse is interrupted (i.e., by jumping to 
dquench = 0) during the period of number loss. Collapse 
processes that were allowed to evolve to completion (until 
N — N remnant ) were not observed to emit jets. Examples 
are shown in Fig. |B| for different T evo ive f° r the conditions 
of Fig. gb. The jet size and shape varied from image to 
image even when all conditions were unchanged, and as 
many as three jets were sometimes observed to be emitted 
from the collapse of a single condensate. The jets were 
also not always symmetric about the condensate axis. 

We believe that these jets are manifestations of lo- 
cal "spikes" in the condensate density that form during 
the collapse and expand when the balance of forces is 
changed by quenching the collapse. We can estimate the 
size of the spikes using the uncertainty principle. After 
a jump to a quenc h = 0, the kinetic energy of the atoms 
in the resulting jet is equal to the confinement energy 
that the spike had prior to quenching the collapse, i.e., 
1-rnv 2 = . •> , where a is the width of the spike in the 
wave function. The anisotropy of the jets indicates that 
the spikes from which they originated were also highly 
anisotropic, being narrower in the radial direction. From 
the widths and the number of atoms in the jets, we can 
estimate the density in the spikes. Plots of the number 
of jet atoms and the inferred density in the spikes versus 
Tevoive are presented in Fig. ^. The jets exhibited vari- 
ability in energy and number that was larger than the 
— 10% measurement noise. 

Overview of the cuegent. . theor etical models. 
Several theoretical papers^lljO'taE^EII have considered 
the problem of collapse of a BEC when the number of 
atoms exceeds the critical number. These treatments all 
use a mean-field approach and describe the condensate 
dynamics usi ng , the . Gross-Pitaevskii (GP) equation. In 
most casest3't2lll§ll3, the loss mechanism is three body 
recombination, but Duine and StooO propose the loss 
arises from a new elastic scattering process. In both 
cases the loss is density dependent and so the loss rate is 
quite sensitive to the dynamics of the shape of the con- 
densate. Since a full three-dimensional anisotropic time- 
dependent solution to the GP equation is very difficult, 
these calculations have used various approximations to 
calculate the time evolution of the_condensate shape. Ka- 
gan, Muryshev, and Shlyapnikovllj numerically integrate 
the GP equation for the case of an isotropic trap and 
large values of the three body recombination coefficient. 
In this regime, the condensate smoothly contracts in a 
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single, collective collapse. Saito and Uedall§ll3 perform a 
similar numerical solution for the isotropic case but with 
smaller values of K3 and observe localized spikes to form 
in the wave function during collapse. Duine and StooftZI 
model the dynamics for the anisotropic case, but use a 
Gaussian approximation rather than an exact numerical 
solution. 

These calculations have all been carried out over a cer- 
tain range of parameters, but none have been done for 
the specific range of parameters that correspond to our 
experimental situations. None of the predictions in these 
papers match our measurements except for the general 
feature that atoms are lost from the condensate. Also 
we see several phenomena that are not discussed in these 
papers. Whether this lack of agreement is due to the fact 
that these calculations do not scale to our experimental 
situation or do not contain the proper physics remains to 
be seen. 

Theoretical Challenges. Collapsing 85 Rb conden- 
sates present a rather simple system with quite dramatic 
behavior. This behavior might provide a rigorous test of 
mean-field theory when it is applied to our experimental 
conditions. Some of our particularly puzzling results are: 

• The decay constant Tdecay is independent of both 



N and \a coUapse \ for |a co H apse | < 100 a , and 
only weakly depends on these quantities for larger 

I ^collapse I • 

• The burst energy per atom dramatically increases 
with initial condensate number. 

• The number of burst atoms is constant versus 

^collapse • 

• The number of cold remnant BEC atoms surviv- 
ing the violent collapse varies between much less 
and much more than N cr , depending on Nq and 
^collapse, but the fraction of remnant atoms, burst 
atoms, and missing atoms are independent of Nq. 

Outlook. From the experimental point of view there 
remain questions to be answered. For one, is the burst 
coherent? It may be possible to answer this by generating 
a sequence of "half bursts" and see if they interfere. For 
another, where do the missing atoms go? If molecules 
and/or relatively high-energy atoms are being created, 
can we detect them? 

It is clear that adjustable interactions opens up a fas- 
cinating new regime for BEC studies. 
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FIG. 1: A example of a ramp applied to the scattering length, 
and a plot of the condensate number versus time after a jump 
to a negative scattering length, a. A typical a(t) sequence. 
ao = 0.529 A is the Bohr radius. The scattering length is 
jumped at t — in 0.1 ms from ai„u to a co u apse , where the 
BEC evolves for a time T evo i ve . The field is carefully con- 
trolled so that magnetic-field noise translates into fluctuations 
in (Lcoiiapse on the order of ~ 0.1 ao in magnitude. The col- 
lapse is then interrupted with a jump to a quenc h, and the field 
is ramped in 5 ms to a large positive scattering length which 
makes the BEC expand. After 7.5 ms of additional expansion, 
the trap is turned off in 0.1 ms and 1.8 ms later the density 
distribution is probed using destructive absorption imaging 
with a 40 /is laser pulse (indicated by the vertical bar). The 
increase in a from a co u aV se to a expan d is far too rapid to allow 
for the BEC to expand adiabatically. On the contrary, the 
smaller the BEC before expansion, the larger the cloud at the 
moment of imaging. Thus we can readily infer the relative size 
of the bulk of the BEC just prior to the jump to a qU ench- The 
density of the expanded BEC is spnlow that the rapid tran- 
sit of the Feshbach resonance poleEj during the trap turn-off 
and the subsequent time spent at B — (a = —400 ao) both 
have a negligible effect, b. The number of atoms remaining 
in the BEC versus T evo ive at a co u a p se = —30 ao. We observed 
a delayed and abrupt onset of loss. The solid line is a fit to 
an exponential with a best-fit value of t co Hapse — 3.7(5) ms 
for the delay. 



FIG. 2: The collapse time versus a co u apse for 6000-atom 
condensates. The vertical line indicates a cr for No = 6, 000. 
The data were acquired with a,mit = a quenc h = (to within 
~ 2 a ). 



FIG. 3: A burst focus, a. Conceptual illustration of a 
radial burst focus, b. An image of a radial burst focus taken 
33.5 ms after a jump from ai n it — to -30 ao for Ao = 15, 000. 
Trad/2 = 28.6 ms, which indicates that the burst occurred 
4.9(5) ms after the jump. The axial energy distribution for 
this burst corresponded to an effective temperature of 62 nK. 
The image is 60 x 310 fim. c. Radially averaged cross section 
of b with a Gaussian fit to the burst energy distribution. The 
central 100 fim were excluded from the fit to avoid distortion 
in the fit due to the condensate remnant (a = 9/xm) and the 
thermal cloud (a = 17/im). The latter is present in the pre- 
collapse sample due to the finite temperature and appears to 
be unaffected by the collapse. The dashed line indicates the fit 
to this initial thermal component. Note the offset between the 
centers of the burst and the remnant. This offset varies from 
shot to shot by an amount comparable to the offset shown. 
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FIG. 4: Burst energies and energy anisotropies. a. and 
b. The axial and radial burst energies versus a co u apse for 
iVo ~ 6, 000 and No ~ 15, 000, respectively. On average, ten 
burst focuses were measured for each trap direction at each 
value of a co iiapse studied. The energies were higher for the 
larger-iVo condensates over the full range of a co u apse studied. 
The vertical and horizontal error bars indicate the standard 
error of the measurements and the uncertainty in a co u apse 
arising from the magnetic-field calibration, respectively. For 
several of the points, the uncertainties are smaller than the 
symbol size and the error bars are not visible, c. The ratio 
of the radial to the axial energies, which is a measure of the 
burst anisotropy. For values of |a co ;; a pse| just past a cr , the 
burst was isotropic for both No — 6, 000 and No — 15, 000. 
At larger values of \a co u apse \, larger- No condensates gave rise 
to stronger anisotropies. For No = 6000, A was 1.6, and for 
No = 15,000, A was 2. When instead we started at ai„u 
■ I OH ao, the BEC was initially more anisotropic (A = 2.4), 
but the burst became more isotropic, with E ax going up by 
-40% and E rad dropping by ~60% for No = 15,000 and 

Icollapse = —100 Oq. 



FIG. 5: Jet images for a series of r evo i ve values for the con- 
ditions of Fig. nib. The evolution times were 2, 3, 4, 6, 
8, and 10 ms (from a to f). Each image is 150 x 255 fj,m. 
The bar indicates the optical depth scale. An expansion to 
o-expand = +250 ao was applied, so the jets are longer than 
for the quantitative measurements explained in the text. The 
jets were longest (i.e., most energetic) and contained the most 
atoms at values of r evo i ve for which the slope of the loss 
curve (Fig. [j]b) was greatest. A tiny jet is barely visible 
for T evo i ve ~ 2 ms (image a), which is 1.7 ms before t co u a pse- 
The images also show how the number of condensate atoms 
decreases with time. The time from the application of a quenc h 
until the acquisition of the images was fixed at 5.2 ms. 



FIG. 6: Quantitative jet measurements, a. The number 
of atoms in the jets versus r evo i ve for the conditions of Fig. 
[TJs. b. The spike density inferred from the kinetic energies of 
the jets. The bars indicate the range of shot-to-shot variabil- 
ity. For the analysis, we assumed the jets were disk shaped 
since the magnetic trap is axially symmetric. The images 
were taken perpendicular to the axial trap axis, viewing the 
disks edge-on. The jets expanded with v ~ 1 mm/s, which 
corresponds to a kinetic energy of ~6 nK and a radial pre- 
quench Gaussian rms width of ~ 0.5 /im. Since the axial 
size was below our resolution limit, we could not measure the 
axial expansion rate. For estimating the spike density, we as- 
sumed an axial width equal to the harmonic oscillator length. 
The atom density in the spikes decreased for larger values of 
|o C oiiop»e|, and was half as large for a coUapse = -100 a as for 
—30 ao- 
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